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ABSTRACT 


A method 1s described whereby a transevaporator 1& used for sampling 
60-100 yl of aqueous sample. Volatiles are stripped from the sample 
either by a stream of helium and collection on a porous polymer, Tenax, 
or by 0.8 ml of 2-chloropropane and collected on glass beads. The 
volatiles are thermally desorbed Into a precolumn which is connected to 
a capillary gas chromatographic column for analysis. The technique is 
shov/n to be reproducible and suitable for determining chromatographic 
profiles for a wide variety of sample types. 

Using a transevaporator sampling technique, the volatile profiles 
from 70 yl of serum were obtained by capillary column gas chromatography. 
The complex chromatograms were interpreted by a combination of manual and 

j 

computer techniques and a two peak ratio method devised for the 
classification of normal and virus infected sera. Using the K-Nearest 
Neighbor approach, 85.7^ of the unknown samples were classified correctly. 
Some preliminary results indicate the possible use of the method for the 
assessment of virus susceptibility. 
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I. INTRODUCTION 


Profile analysis has been the subject of considerable investigation, 

with various degrees of success, in food quality evaluation, in 

air pollution studies, in the characterization of automobile 

exhaust gases, marine sediments, engine oils,^^^^ and marine 

oil spills. The object of gas chromatographic profile analysis 

is to correlate a characteristic "finger print" pattern in a chromatogram 

with the qualitative and quantitative properties of a sample. 

In recent years, this approach has become increasingly important in 

the biomedical field for the diagnosis of disease and the study of 

metabolic disorders. A complete profile of all the constituents 

of a biological fluid is at present impossible with available analytical 

techniques. The goal at which most workers aim is the development of a 

(30 

complete profile of a selected group of substances (eg. organic acids, v * 

31,36,44,45) acids, Steroids, carbohydrates, ^ 

(19-21 23 24 2fl 37 

etc.) or of compounds with similar properties (eg. vola, tiles) ">>>-> 

A comparison is then made between the profile obtained with "norriial 
samples" and "pathological samples" to establish any quantitative 
differences that might be of value for diagnostic purposes. A knowledge 
of the chemical constitution of abnormal peaks or peaks of abnormal 
concentration may suggest a possible biochemical reason for the disease 
and lead to new methods of treatment. 

Except in a few cases in which a strong clue exists to implicate the 
role of causative agents to a disease disorder most work still relies on 
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a serendipitious finding in a screening program involving the study of 
many diseases. To maximize the chance of discovery as many substances 
as possible are included in the profile. Capillary columns are generally 
used for the analyses due to their superior resolving power. 

Volatile metabolites excreted in urine have been extensively 
studied by Zlatkis Liebich et al >37,51 ,52) 

Matsumoto and Novotny . Abnormal excretion of C^-Cg alcohols, 
octanols, cyclohexanone, 2-heptanone and pyrazines occur in cases of 
metabolic disorders of diabetes mellitus. 

Clinical laboratories depend to a great extent on blood as opposed 

to urine for diagnostic purposes. Little data is available concerning 

the occurrence and concentration of volatile substances in blood. 

Technical problems such as the availability of sufficiently large samples 

and the protein content of blood serum have rendered their qualitative 

T53l 

and quantitative determination more difficult. Dowty et / used 
Tenax GC to analyze volatiles in plasma of patients 'with chronic renal 
failure. Stoner extracted volatile substances from 15 ml of 

plasma with ether. Liebich and Woll^^^^ analyzed 5 ml of serum by GC-MS 
and mass fragmentography using a headspace sampling technique. In all 
these studies the sample size used for analysis ranged between 5 ml and 
55 ml and the chromatograms produced contained relatively few volatile 
components compared to urine. 

The components of interest are generally present in low concentration 
with the result that special techniques are necessary to overcome the 
problems associated with the small sample size, low sample capacity of 
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the capillary column and the need to introduce the sample as a con- 
centrated plug. Two techniques that have been widely used for the 
concentration of volatile substances from dilute aqueous solution are 
solvent extraction^^’^’^^ ^ and headspace analysis 
using a suitable adsorbent. The 

techniques are complementary for samples that contain organics covering 
a wide volatility range. A shortcoming of the' extraction technique is 
that some early peaks are masked by solvent whereas poor representation 
of late peaks occurs with the headspace technique. The extraction may 
be performed using a micro-technique but the headspace method generally 
yields unsatisfactory results where only small volumes of sample are 
available. 


Zlatkis and coworkers introduced a direct micro-extraction method 
for profile analysis using only 25-100 yl of serum. Zlatkis and 
Kim have recently described an improved method of micro-extraction 
using a transevaporator. An aqueous sample is placed on a suitable 
adsorbent packed in a small tube, the organics of interest are extracted 
v/ith a suitable solvent and transferred to a collection tube packed with 
glass wool and excess solvent is evaporated. The volatiles are then 
transferred to a capillary column via thermal desorption into a precolumn. 
The extraction procedure takes less than 30 minutes and may be used with 
as little as 25 yl of sample. The transevaporator extraction is about 
three times more efficient than a direct extraction. Purified 2-chloro- 
propane was found to be a satisfactory extracting solvent and Po’rasil E 
was found to be the best adsorbent out of 6 possible candidates. 
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The volatiles are present at trace levels in biological fluids v;hich 
consist principally of water so that the sampling technique has to serve 
as a concentration device capable of reproducibly stripping a diverse 
range of compounds from a v/ater matrix. The original "transevaporator" 
method has several shortcomings. A small excess of moisture sometimes 
causes plugging of the analytical capillary column, and the chromatogram 
was lacking in headspace component.,;. The recovery and reproducibili ty 
of standard compounds were also lower than expected. It was necessary 
to modify and improve this "transevaporator" sampling system for its 
general application to biological fluids. Good precision was required 
and a more general profile containing "headspace" and "solvent 
extraction" information from the same sample. 

When one is dealing with the metabolic disorders that lead to gross 
metabolic changes, there is no need for computer aided evaluation. 

However some clinical disorders may result in more subtle, yet perhaps 
important, changes in the metabolic pattern. The very nature of this 
approach leads to problems of an analytical nature. For the case of the 
analysis of volatiles, the current sampling techniques generate very 
complex chromatograms which only reveal a full picture of their complexity 
when high resolving capillary columns are used for the gas chromatographic 
separation. The amount of information subsequently generated can no 
longer be handled by empirical means and the use of pattern recognition 
techniques and computer sorting are obligatory. ^ Although the 
selected sample can be characterized as "volatile" under that heading 
is masked the chemical complexity of the mixture which covers the 
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complete spectrum of polarity. The physical property shared by all 

components is that they fit into a distinct boiling point range. The 

successful use of such techniques as pattern recognition makes high 

demands on the reproducibility of the profile; a feature which it is all 

too easy to demonstrate is more often impaired by poor sampling techniques 

( 10 ) 

than chromatographic error.' ' 

Various pattern recognition techniques have been described to 
correlate gas chromatographic profiles with the properties of a sample. 

The purpose of pattern recognition is generally to categorize observed 
experimental data as being characteristic of a particular sample type. 
Moskowitz^^^ used multi linear regression analysis to yield an overall 
dissimilarity of odor quality difference in flavor assessment. Ours 
^ used a linear learning machine pattern classifier to analyze 

low resolution mass spectra. The linear learning machine performs as a 
binary pattern classifier - yes or no, plus or minus, or class I or 
class II. A classification vector is computed by an interactive 
training data set using known spectra that may overlap in' any or all 
components. This technique has a disadvantage for the non-Tinear multi- 
pattern data set compared to the K-Nearest Neighbor technique (KNN). 

The K-NN classification method^^^"^^^ is a multi-dimensional classifier 


with no increase in complexity and is computationally and conceptionally 
quite simple. An unknown pattern is classified according to the 
majority vote of its K-Nearest Neighbor (K = 1 ,3,5. . .etc, ) in n-space. 
Computationally, all that is necessary is to calculate and scan a 
distance matrix between an unknown point and points of a known group 
established from the analysis of a training set. 
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Profile analysis of volatile compounds in body fluids could be a 
significant part of a large-scale preventive medicine program provided 
that an accurate and long-term reproducible analytical system could be 
developed. Specimen collection, storage and transportation should be 
controlled carefully in order to eliminate any artifacts. The combination 
of high resolution capillary column gas chromatography and computer aided 
pattern classification systems, opens up new possibilities not only for 
the study of gross alterations but also for the detection of more subtle 
changes in profiles. 

In an effort to establish the methodology for the early diagnosis 
of viral diseases; a series of capillary column chroma togrpahic profiles 
were obtained using a modified "transevaporator" sampling technique. 

The complex chromatographic profiles were analyzed by computer techniques 
to develop a means of differentiation between normal and virus infected 
sera using a two peak ratio method. A training set was developed and 
the K-Nearest Neighbor technique used to establish the predictability of 
virus infected serum identification. Gas chromatography-mass spectro- 
metry was used to identify the inost prominent serum components important 
to this study. 
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II- EXPERIMENTAL 
A. Adsorbents and Reagents 

For trapping and concentrating the volatiles from biological fluids, 
the adsorbent should provide adequate selectivity, thermal stability and 
surface activity. 

Porasil E (80-100 mesh) and glass beads (80-100 mesh) viere obtained 
from Analabs, North Haven, Connecticut. Porasil E is in the form of 

perfectly spherical beads of porous silica with a specific surface area 

2 f73l 

of 25 m /g. According to Cadogan and Sawyer,'' ' upon heating from 

200°C to 500°C silica loses all hydrogen-bonded water which provides the 

active sites for adsorption. Since Porasil E was conditioned at 300®C 

prior to use, water in serum was retained by Porasil E very effectively 

after sample elution through this small column. 

Tenax GC (60-80 mesh), a porous polymer of 2, 6-di phenyl -p-phenylene 
oxide, was obtained from Applied Science, State College, Pennsylvania. 
Tenax GC has been characterized as having unusual and desirable 
chromatographic properties. It has been shown to fulfill the 
requirements of efficient adsorption and thermal desorption for higher 
molecular weight compounds. The temperature stability of Tenax GC is 
relatively high (up to 400°G), and it is widely used for headspace 
sample analysis. 

Glass wool (Pyrex brand) was obtained from Corning Glass Works, 


Corning, New York. 



8 


All of the above were washed with distilled ether and conditioned 
at 2S0°C in a stream of dry helium for 12 hours. 

2-chloropropane (Eastman-Kodak, Rochester, New York) and n-pentane 
(J. T. Baker, Phillipsburg, New Jersey) were distilled from phosphorous 
pentoxide (J. T. Baker, Phillipsburg, New Jersey). Diethyl ether 
anhydrous (Mallinckrodt, St. Louis, Missouri) was distilled from lithium 
aluminum .hydride (Alfa Products, Danvers, Massachusetts) and water was 
distilled from potassium permanganate (Fisher Scientific, Pittsburgh, 
Pennsylvania). 

B. Serum Samples 

A total of thirty-six serum samples from twelve male volunteers who 
had been exposed to either "England" or "Rhinovirus" were supplied by 
Dr. H. J. Schneider (NASA, Houston, Texas, U.S.A.). These samples 
form part of a study of immune response to respiratory virus infection 
conducted by Dr. R. Couch, Dr. V. Knight and Dr. S. Criswell of the 
Influenza Research Center (Houston, Texas). Serum samples (200-800 yl) 
v^ere stored in glass vials at -20°C prior to analysis. 

The serum samples fall into two categories. The first (Nos. 1 ->-6) 
were obtained from volunteers who developed clinical symptoms of influenza 
after virus infection. The second group (Nos. 7 ^ 12) were obtained 
from volunteers who did not develop clinical symptoms after virus 
infection. Each volunteer provided three serum samples corresponding to 
a baseline sample taken prior to virus infection (Group I), a serum 
sample taken one day after exposure to the virus (Group II) and a final 
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sample taken 14-21 days after virus infection (Group III). Group (I) 
represents normal serum samples. Group (II) virus infected serum samples 
(irrespective of whether clinical symptoms developed after exposure) and 
Group (III) virus infected serum samples from which all volunteers had 
recovered clinically from the infection. 

C. Transevaporator Sampling Hethod 

1, Sampling Apparatus 

The transevaporator used is shown in Figure 1 and is essentially of 
the same -design described by Zlatkis and Kim.^^^^ Significant differences 
are that the dimension of the Porasil E tube is changed to smaller size 
(70 mm X 2 mm I.D.) and that a water condenser is included. 

The micro elution column was a glass tube containing 0.3 ml of 
Porasil E with a plug of glass wool at both ends. Prior to packing, 
Porasil E v/as conditioned again at 300®C by purging with nitrogen at 
20 ml/min for 3 min. The glass bead trap (3.4 ml of glass beads) and 
Tenax trap (1.8 ml of Tenax GC) serve as a concentrating ‘column, in a 
glass tube (110 mm x 10 mm O.D., 8 mm I.D.), tightly packed with glass 
v;ool. The elution apparatus consists basically of a custom made 
"transevaporator" and a conical tip centrifuge tube with a glass stopper 
(§ 10/18) as shown in Figure 1. “ 

2. Sampling Procedure 

Figure 2 is a schematic diagram of a combined headspace and 
extraction sampling procedure. The sampling procedure is performed in 
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two steps, the first being collection of low molecular weight volatiles 
on a Tenax tube, the second being the solvent extraction of volatiles 
coupled with collection on glass beads. 

In the first step the sample (25 to 200 pi) is introduced to the 
inner tube which contains 0.3 ml of Porasil E, a strongly hydrophilic 
porous silica adsorbent. The upper tube contains 1.8 ml of Tenax, a 
2,6-diphenyl-p-phenylene oxide porous polymer, which is an excellent 
trapping material for headspace analysis. A condenser is included, in 
the apparatus to prevent excess water vapor from reaching the Tenax trap. 

A stream of helium gas (16 ml/min) is passed through the apparatus for 
5 to 10 minutes to transfer the volatile compounds from the porasil 
tube to the Tenax tube which is then removed from the apparatus. Helium 
(16 ml/n,1n) is passed through the tube for about 3 minutes to remove 
traces of water which v/ould block the precolumn during desorption.: There 

is no apparent change in the resultant profile due to this procedure. 

The Tenax tube is removed from the apparatus and placed in a desorber at 
280°C for 10 minutes. A stream of helium (7 ml/min.) desorbs the volatiles 
into a precolumn, cooled in liquid nitrogen, which is then connected to 
an analytical capillary column for chromatography. 

In the second step the apparatus is converted to the form described : 
by Zlatkis and Kim^^^^ by removing the condenser. The Tenax tube is 
exchanged for a tube of identical dimensions containing glass beads. 

0.8 ml of 2-chlorDpropane is introduced into the tube and is swept through 
the porasil tube into the upper tube by applying gas pressure. Soluble 
organics are extracted and transferred to the glass beads while most of 
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the water and high molecular weight material is retained by the Porasil 
E. The Porasil E tube is removed, the apparatus is reassembled and the 
outer tube is placed in an air bath maintained at SO^C to heat incoming 
helium which is passed through the apparatus at a rate of 16 ml/min for 
10 minutes to remove excess 2-chloropropane and traces of water. As in 
the case with the Tenax tube there is no apparent change in the resultant 
profile due to this procedure. The organics are transferred from the glass 
bead trap to the capillary column by the same technique as is used for 
the Tenax tube. 

3. Transference of Sample from Trapping Tube to Analytical Column 

The trapping tube was placed in a desorption chamber maintained at 
280°C, connected to a stainless steel precolumn (30 cm x 1 mm I.D. , coated 
with SF-96) and cooled in liquid nitrogen. The volatiles were transferred 
to the precolumn by passing helium through the system at 7 ml/min for 10 
minutes. The precolumn was then connected to the analytical capillary 
column, the first 30 cm of which was cooled in liquid nitrogen. The 
sample was transferred from the precolumn to the analytical column by 
heating the former to 180°C with an air heating gun for 1 minute while 
passing helium through the system at a rate of 1.5 ml/min. 

D. Capillary Column Gas Chromatography 

1. Analytical Capillary Column 

In this study, all analyses were carried out on stainless steel 
capillary columns (100 m x 0.25 mm, I.D.) which were purchased from 
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Handy and Harman Tube Co., Norristown, Pennsylvania. Newly acquired 
columns were first checked for leaks and then modified at both ends by 
slipping short pieces of stainless steel tubing (40 mm x 1.0 mm I.D. and 
1.6 mm O.D. ) over the end, after silver soldering, the sleeve was fitted 
with 1/16 in. ferrule and nut. 

Cleaning of the new columns was done by washing with several 25 ml 
portions of acetone, n-pentane , chloroform, alcohol, water and acetone 
at 500 psi . After drying the column at 150®C for one hour, the quality 
of the empty column was tested by injecting a series of test substances 
(cyclohexane, benzene, diethyl ketone and ethanol) individually at SO^'C. 
All test substances should be free of significant tailing. To improve the 
quality of empty columns which show, severe tailing, the column was filled 
with etching solution (5 g ferric chloride, 15 ml cone. HCl and 60 ml 
distilled water) for 2 hrs and washed -with distilled water several times, 
followed by solvent washing. 

The liquid phase used in this study was 5% Witconol LA-23 (Witco, 

Houston, Texas), an alkyl aryl polyether alcohol, in CH 2 CI 2 solution 

% 

(W/W). A post column of identical I.D. is attached at one end of the 
analytical column and the other end is immersed in a beaker filled with 
water. The flow rate throughout the coating process was controlled at i 
one bubble per second. After drying overnight at room temperature by 
passing nitrogen carrier gas the coated column was conditioned in the 
gas chromatograph at a low flow rate and temperature programmed at a 
rate of 0*5®C/min from 40° C to 160°C. 


15 


The efficiency of a coated column v/as checked with a decalin mixture 
(decalin:diamyl ether =1:1, V/V) at 120®C, isothermally. The resolving 
power was also tested with an aromatic mixture (cyclohexane, benzene, 
ethylbenzene, meta- and para-xylene) at 80®C, isothermally. Most of the 
columns which were coated with this method gave 250,000 theoretical plates 
and provided base line separation between meta- and para-xylene. 

2. Gas Chromatographic C onditions 

A Hewlett-Packard 5830 A gas chromatograph (Hev/lett- Packard, 

Avondale, Pennsylvania) equipped with flame ionization detectors and 
connected to a Hewlett-Packard -18850 gas chromatography terminal was 
used. The serum volatiles were separated on a stainless steel (100 m 
X 0.25 mm I.D.) capillary column coated v/ith Witconol LA-23. With a 
helium carrier gas flow rate of 1.5 ml min \ the column v/as held iso- 
thermally at 50®C for ten minutes and then programmed at 1.5°C min""* 
to 160°C and maintained at this temperature for eighty minutes. 

3. Gas Chromatography-Mass Spectrometry (GC-MS) 

Gas chromatography-mass spectrometry (GC-MS) was performed on an 
LKB 9000 instrument (LKB-Prodkter AB, S-161, 25 Bromma 1, Sweden) with 
a single jet separator and a Perkin-Elmer 900 gas chromatograph. The 
ionization voltage was 70 eV and the scan time for a mass range of 20 to 
300 atomic mass units was 4 seconds. The temperature of the separator 
and ionic source were set 220°C and 250°C respectively. Analytical 
procedures and chromatographic conditions were the same as described in 
the previous section. When possible identification was confirmed by 
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comparison with standard compounds available in the laboratory, otherwise 
manual interpretation and comparison to library spectra^^^*^^^ v;as used. 

E. Data Handling and Computer Interpretation 

The operations which constitute the interpretative procedure are 
summarized in Figure 3. Visual inspection was used to identify those 
peaks common to all chromatograms. The peaks were normalized and this 
formed the data set for calculation of system variations and variation 
between individuals within each group. A second data set selected on 
the basis of the magnitude of their difference between groups was used 
to calculate all possible peak ratio combinations and diagnostic ratios 
selected. These ratios were autoscaled and used as input for K-nearest 
neighbor classification 

1. Autoscalinq^ ^^*^^^ 

In autoscaling, each dimension was adjusted to a mean of zero with 
unit standard deviation. Hence: 

y ^ _ T m m 

im - 0 ^ 

where is the mean and is the standard deviation of the value 
for the ith feature. is the new autoscaled value. This preprocessing 
minimizes the effect of differing data magnitude. It was known that 
autoscaled data gives high rates of correct classification in’K-NH method.- 

2. K-Nearest-Neiqhbor (KNN) Classification ^^^*^^"^^^ 

This pattern classification technique is based upon the assumption 
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that "the closer together" patterns are, the more closely related are the 
sets of data. Closeness was defined as the distance between an unknown 
data point and points of a known group (training set). The distance 
between point 1 and m is defined in this case as a Euclidian distance: 

^Im " = 1 . 

where k is the number of dimensions in factor space. In the case of the 
single (k = 1) nearest neighbor (1-NN), the unknown point is classified 
as the group which gives the smallest distance from the points used for 
the training set. 

3. Computer Programs for Profile Analysis 

All programs were written in Fortran IV and run on the University 
of Houston computing center's Honeywell 66/60. These computer programs 
consisted of a short main program which inputs a set of data, initializes 
a number of parameters, and calls the several subroutine subprograms; 
percent normalization, calculation of relative standard deviation, 
writing, calculation of all possible combinations of two peak ratios, 
autoscaling, sorting in ascending order and classification by K-nearest 
neighbor methods. 

The array A contains the raw integrated peak area. It is dimensioned 
for up to 36 chromatograms (N) of up to 37 selected peak areas (H) per 
chromatogram. The following is an example of an actual computer program 
used in this study. 


on o c o o o n n n o n o r» {-> o o o o o c*> o n n n n 
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****** A *********** **/i'*****A*»t****AA**** A ****** A********!*; 

AUG. 10 » 1978 K.y.l.EF- 


PROFILF. ANALYSIS OF BIOLOGICAL FLUID 

N « NdiMPER UF COLUm.-nIS 'iS^-MPLE nUMBE’RS) 

M = NUOBER OF WQ'.'-S (SELECTED PEAKS OR PEAK RATIO) 

A = Data points ARRAY: A(M,N) 

S = PERCENT normalized DATA POINTS; SCM.N) 

R = PEAK RATIO; RCM,N) 

4SLB = AUTOSCALED POINTS ; ASLBCM.N). 

DOliPLE PRECISION 8 , a SLB , meanB , SOQ , PTCV ,R 
DIMENSION A(3V,36) ,3(37,56) ,ASLBC37,36) ,meaNBC37) , 
1 SD8C37) ,PTCVC37) ,R(7,36) 

READCS,100) N,M 
initialise DATA POINTS- 

DO 10 J = I , N 
A(1 ,J)=0.0 
&(I,J)=0.0 
RfI,J)=O.0 
ASLb'Cl , J ) :::0,0 
MEANS(i)=0,0 
SD8Cn:^0,0 
PTCV(x)=0.0 
10 CONTINUE 

'■ WRITE (6,118) 
write. (6,1-lv) 

DO 15 I = 1 ,M 

read (5,101) CACI,J), J=1,N) 

WRITE C6,i2l) I, CA(I,K), K=1,N) 

IS CONTINUE 

NORMALIZATIONCPERCENT) OF DATA 
C.ALL NORM (A,N,M4B) 

CALL WRITE (B,N,M) 

C CALCULATION OF REPRCDUCIBIUI TY 
WRITE (6,330) 

L = l 

NN 2 12 
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16 , 

KL 3 

CALL CALC CB ,L,N,NN«KL,M,M£ANB,SCB ,F'TCV ) 

WRITE (6,340) 

L - 2 

call calc C8,L,N,'NL',KL,M,m£anB,S0B,PTCV ) 

WRITE C6,3B0) 

L = 3 

CALL CALC CB,L,N,L'N,KL,M,MEAMB,S08,PTCA’ 5 

100 FORMAT (2110) 

101 FORMAT ( 3(VF8,f),/) .9F8.0 ) 

113 format ■•**•*« DATA PO.T.N'B ************* 1 !? 


l?V*yf*******v#i-**>>**il<********'3r*iir'*'*'rtA , .'/) 


119 format 

(// 
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1 "10-3", 3X, 

1 "U“i" ,3X, "11-2" ,3X," 11-3% 3X,» 12-1", 3X, 

1 "12-2" ,3X, "12-3" , /// ) 

121 FORMAT (13, 13F7.0,/,3X, IfiFT.O, // ) 

330 FGRMATCIHI ,//,'%*BHFORE INFECTION 12 SAMPLE S** ",// 5 
3^10 FORMAK Irtl , //, "**aFT£R INFECTION 12 SAMPLES* A-" ,// ) 
350 FORMATClHl ,/./,’%*AFT£R CURE , 12 SAMPLES**'',//) 

STOP 
END 
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Si- * * 5^ * * * * + •<;• ■* * iJr A- * A A * * V * < )V ■#>■ * A '^ * * * * ir A * A -A V: * ^ * y- 4 * •*• ^ * <r * 

•AAAA'AA>’r>l:'A**A><i<,"AA'«>'cA'A*+AA-A*AA*AA'A***Ay:*<fk'*'AA:<*'A*'i»*'fev>lt** 

pnoLED VI PUS sample, plprooucibili rv 

NOHMALlZATIOfM , 6 SAMPLES 
PEAK PATIO, ANQ MEA^J . SO, PT, COEF. VAR. 

kA'A*'*>'kkk*k>f*'icT<rk^*kk'>it*A*AkiSfk**#rAkkkk*****k*kiir**x^k*kkk* 
**?!:*»i(*<r**kk*A **k*kk**A Akk***^!** krtkkk***** *'^*'4!**A«kkjir*k'*r 


AUG. 1 0 , 1978 K.Y.LEF, 


PROFILE analysis OF BIOLOGICAL FLUID' 

N r. NUMBER OF COLUMNS (SAMPLE NUMBERS) 

,M = NUMBER OF RuivS (SELECTED PEAKS) 

A s data points array; a(m,n) 

B - PERCENT NORMALIZED DATA POINTS; B(M,N) 

.ASLB - AUTCSCAlF.D points ; ASLB(M,N) 

R = PEAK RATIO 

DOUBLE PRECISION B , ASLR ,M£AWa ,SD6 ,PTCV ,R 
DIMENSION. A(37,7) .B(37,7) , A SLR ( 37 , 7 ) , ME A NB ( 37 ) , 
1 SD6(37) ,PICV(37) ,R(a,7) 

C 

REAOC5,100j N,M 
C INTTT.au ZE DATA PCIMTS 
00 1 C I = 1 , M 
00 10 J = 1 , N 
A(IU3=0,0 
8( I , J}=0 . 0 
H(1 , J) =0.0 
ASLBd ,J)=0.0 
MEANB(I)=0.0 
SDB(n=0,0 

PTCV(I)=0.0 •• d . 

10 CONTINUE 

c 

WRITE (6, IPO) 

DO 20 1= 4-, M 

read (5, ion (A(I,J)t J=1,N) 

WRITE (6,121) I , (ACUK: ,K = 1 ,N) 

20 CONTINUE 
C 

C NGRMAL!2ATIQN(PERCENr) OF DATA 
CALt. NORM ( A,N,M,8 ) 

C CALCULATION Or MEAN, S.D., P£«Ct NT COEFF, VAR. 

L = 1 

MM s 6 
KL=1 

CALL CALC (fi,L,N,NN,KL,M,M£ANR,SDB,P7CV ) 

C calculation of peak ratio Or POOLED SNSUM 
WRITE (6,122) 



22 


16.^126 



DO 

23 

1 

= 1 
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DO 

25 

J 

= 1 

, N 





IF 
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u , 

EQ. 

6) 

GO 

TO 25: 


IF 
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) GO 

TO 

28 


IF 
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) GO 

TO 

29 
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GO 

TO 

25 






2R 

R(I 

t 
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CRCI 

, L) « 


23 CONTINUE 

C CALCULATION OF ME AN S . D , , PEPCEN'T COEPr. VAS, OF RATIO 
M-2 

CALL CALCI (R,L,N,NN, KL,M,M£AMB, SOB, PTCV ) 

C AUTOSCAI.ING OF PEAK RATIO 

CALL AUTS CR,N,M,MEANB,sDB,ASLB ) 

C calculation Or .MEAN, SO, PT , COEF , V Af< . OP AUTQSCALED POINT 
CALL CALCI CASLB,L,N,NN,KL,M, MEANS, 5Df3»PTCV) 

100 FOPMA.TC2I105 

101 FORMAT ( 7F10‘.0 ) 

120 format < ' I ' ,SX , "**<<0ATA POlMTS**)t'\/, 1 1X,”PEAK MQ” , 

i 5X,"SAMP 1«, 6X, "SAMP a«,8X, ''SAMP 3''*,3X, 

1 "SAMP ax, "SAMP 5",SX,»SAMP 6",8 X,"SamP 7“,. 

n ) 

121 FORMAT (/,lOX,I5,7Fl^,O0 

_122 FORMAT (IHl ,///,aX,"**PEAK RATIO OF i/ii, 5/7 **** 

1 •',//, 1 1 y , "PEAK NO"» 

i "SAMP a", ax, "SAMP 5'‘«ftX,"SAMP 6",6X,"SAMP 7", 

i 5X,"SAMP l",8X,"SAMP 2",BX,"SAMP 3",8X, 

1 // ) 

123 format C 10X,I5, 7Fia,2, // D 
STOP 

ENO ' 



i ri o o r; o r> o o n o <■> n n n o r» n o o o 


i * il- * A * A * * ****** * 5 ;: * * A * V: it; >> * * -A: ♦!<:*■*:*★* * 'If * X ****** Vr * * * * >;* 

***************************** ,S ,; * * i, ;V ************ K * X ***** * 

PROFILE AMALY3ir> tDF hf£F OR F. , A P TE R Ani) CURE VIRUS SAMPLE 
SELECTED 7 PE AK s c FORMAL J Zr. D ) P 16 ::I, P1R = 2, 

P20--3, P25 = a, P29nS, P3J:i6, P:s^ = 7, 

RESPECTIVELY Af^D C ALC lil. A T IO>N OF PEAK RATIO 
REPRODUCIBILITY Or PEAK RATIO WiTHIiM CROP 

>f*****'A'******K ************************************ ****** 
*************************************************** **x** 

AUG. 10 y 1976 K.Y.LEE 

PROFILE ANALYSIS OF BIOLOGICAL FLUID 

N - number of columns (SAMPLE NUMBERS) 

M s NUMBER OF KO‘vS (SELECTED PEAKS) 

B = PERCENT normalized DATA POINTS; B(M,N) 

ASLvS p AUTOSCALED POINTS; ASLBCM^N) 

R = peak ratio ? R(M,N) 

DOUFLE PRECISION 8 , A SL 3 , R , ME A N6 , SD8 , PTC V 
DIMENSION B( 7.36) ,A3LPC7,-U0 ^RCT.Sb) , 

1 M£ANB(7) ,SD6 (7) ,RTCV(7) 

READC5, 100) N,M ’ • 

C ■■ INITIALIZE DATA POINTS 
DO J vO I = UM 
DO i 0 J P i , N ■ 

8(I,J)=0,0 
R(I , J)=0.0 
ASLBCI .J)=O,0 
ME.4NHCI)p0,0 
SDsCDpO.O 

PTCV(I)=0.0 ■ ’ . . 

10 CONTINUE 

c 

WRITE ( 6 , 120 ) 

DO 1'5 I - 1 ,M 

READ (S.iO-i-) (B(I,J) , J=1,N) 

WRITE (6,121) I, (6(I,K), K=1,N) 

15 CONTINUE • 

C 

write ( 6 , 220 ) 

CALL WRITE <3,N,M) 

C CALCULATION (jf PEAK RATID(NORMALIZEO) 

CALL RATG(5 ,N ,M,R) 

100 FORMAT(2I103 

101 FORMAT ( 3(10F6.2,/), 6F5.2 ) 

120 format ( 1 hi , // ,6X , ''*v.seLECTF-D 7 PEAKS, NORMALIZED 

1 ******* ***K*** ***** A* ****************,*****'* y // / , 

1 •‘P£AK’‘,2X,"l-rS^X,"l-2'^,0X,/’l"5‘',ax, 
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16.431 


121 FOR/M T 
220 FOF.MiT 

1 *A**;if * 

STOP 


"2«1",4X,"2"2”,4'^,"2-3",4X,’'3-1",4X."3-2'' , 

4X , "S-S*" , ax , "^-r' , . ''4-2" , 4X , "^-3” , 'IX , 

5 “ 1 " , 4 X , " 5 - 2 " , a X , " 5 - 3 " , 4 X , " 6 - i " , a X , " 6 - 2 , 
4X,"6-3'' ,4X,//,6X,''7-t%4X,''7-2",aX,''7~3", 
ax . ''8~r' , 4X . ''8-2” ,4X , ''8-3" , 4 X, "9-1 " ,4X , 
”9«2’S 4X , "9-3'' , 4>; , " 1 U“1 " ,3X , " 1 0-2" ,3X . 
"10“3" ,3X. 

'* 1 1 - 1 " , 3 X , " 1 1 - 2 " . 3 X , " 1 1 - 3 " . 3 X , " 1 2 - 1 " , 3 X , 
"12-2",3X^»l2-3'' ♦ /// 5 
CI3, 18F7 .2,/,'iX, 16F7.2. ///) ‘ 

ClHl ,//, "*vc**-A^fJRM.aLIZEO OATA, 7 SAMPLES 

**it***:ic* Ik ******************************* ".///) 
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PROFILE ANALYSIS 

AUinSCALIN'G GF SELFCTFD 2 PEAK RATIO, . 1/A AND 5/7 

^(rV.'*^'T»n»r;fc*ii.'*Vr*x*r3if*A**Vr5!:'ff**^*->t*>s'*i!r****yc-*ivi'^*yr**Kf***w****»*** 

lit,V*-**^*<r*Ji(*‘/rjt********-<!.'**:if**V'r*Tif*ilr******-*'*************7ir** 

AUG, 10, 1978 K.Y.LEE 


PROFILE analysis OF BIOLOGICAL FLUID 

N = NUMBER OF COLUMNS (SAMPLE NUMRF-RS) 

M = NUMBER OF RO^S (SELECTED PEAKS OH PEAK RATIO) 

R = PEAK RATIO data; RCM,N) 

ASLB = AUTOSCALED POINTS ; ASLB(M,N) 

DOUBLE PRECISION R , A SLB , ME ANK , SOB , P TC V 
dimension R(2,3S) , ASLBC2,36) , ME ANB ( H ) , SOB (2) , PTC V ( 2 ) 

REAO(5, 100) -N,?^ 

INITIALIZE DATA POINTS 
' DO 1 0 I = 1 
DO 10 J = 1, N 
RCI , J)=0.0 
ASLBd , J)=0.0 
MEANB(I)=0,0 
SDh(I)nO,0 
PTCV(I)=0,0 

10 continue 

WRITE (6,120) 

DO IS I = 1 ,M . 

READ (5,101) (R(I ,J) , J = 1 ,N) 

WRITE (6,121) I, (R(I,K), Kc:l,N) 

15 CONTINUE 

WRITE PEAK _RATIQ OF EACH GROUP 

write C6,1S1) 

CALL WRITE (R,N,M) 

C calculation Or MEAN, S.D,, PERCENT COEPT. VAR, 

L = 1 
NN = 36 
KL = 1 

CALL CALC (R,L,N,NN,KL,m,M£ANO.SDB,?TCV ) 

C AUTGSCALING OF PEAf< RATIO; Mf£AN::0.0, S.D.-l.O ,, 

CALL AUTS (H,N,M,ME«NB,SDri,A3LS ) 

WRITE (6, 126) 

CALL WRITE (ASL3 ,N,M ) 

c ■ 

100 FORWAT(2IIO) 
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101 F0T<MAT ( 5(iCF3,2,/), 6F8.? ) 

120 FOf<N’AT ClHl SEI.FCTFO 2 PEAK RATIO**-”, 

1 "I/« A^iO 5/7, *-****rr^*A^-*^*:VJ^******'' ,//, 
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121 POPf^-AT Cn, 18F7.2,/,3>', 16F7.2, //////) 

■'l26 FOPMAT C 1 ' ,8X , ''**auTQJ?CAUIA>G OF PEAK RATIO*-***** ’ 

1 lir**'*1!r**-->r1< **•;(*•* lV**il;**y****i*-ilr'i:*>l:*ik****-*******-*^l<”,/ / 5 

151 FGPV.AT ni-Jl,//, "** SELECTED 2 PEAK RATIO ********* 

l***'*****?K<:****«r*t(rA*>;*******ii<*W»nir^**A*******i»'***',//) 

STOP 

END ■ - - 
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**iiyt<t***if*-ir*Ai-it:*-4-****'i*-»fii(***ir*Ti+**-Jt*'*r***^j>-#*****^i!f*iMf>* 
* * * * •> * * A- -i *• * * * ■* w * t * * * ■;► * if lif ^ * )iif * * * (if If -if * >» * * * * * A * * *• A * * :V * 


K*’NiN PATTERN Rt'CQGwITlt]'-J 
NORMAL OR VIRUS INFECTED SERUM 


* A * * ■if >t * ' 

AA***AA*Ai 


AAA*******: 


*i-Ai**AA*A*AAAAA*»:* 
r * * A A 


■ AUG, 10, 1R78 K.Y.LEE 

PROFILE analysis OF BIOLOGICAL F LlJ I D K-NN_ CL A SSI F I C A T ION 

M K NO, OF UNKNOWN SAMPLES C ROW) : M(I) 

L = NO. OF training SETCROWI ; LCJD 
'N = NO. OF FACTOR SPACE (COLUMN) ; N(K), K=3 
% = UNKNOWN sample SET; YCI,K) 

Y - KNOWN TRAINING SET ; YCJ,K) 

D = euclidean DISTANCE IN N ( K ) SPACE BETWEEN PDrNTI,J 
double precision X,V,D 

DIMENSION XC26,2), YC10,2), 0 ( 26 , 1 0 ) , JJ C 1 0 D 

C 

READ (5,100) L,M,N 
C INITIALIZE DATA POINTS 
DO 10 I = 1,M 

DO 11 J = 1 ,L 

DO 12 K - 1,N 

XCI,K)=0.0 
Y(J, 10=0,0 
DCI,J)=0,0 
12 CONTINUE 
n CONTINUE 
10 CONTINUE 

c 

WRITE (6,110) ■ ■ 

DO 15 I = 1 ,M 

read (5,105) (X(I,K) , K=i ,N) ■ 

WRITE (6,112) I. (XCI,K), K=1,N) 

15 CONTINUE 

WRITE C6,1T3) 

DO 17 J = 1 ,L 

READ (5,106) CY(J,K), K=1,N) 

WRITE (6,112) J, (Y(J,K), K=1,N) 

17 CONTINUE 
C 

C CALCULATION OF EUCLIDEAN DISTANCE IN' 2 FACTOR SPACE 

WRITE (6,220) 

0 0 2 0 I = I , M 
DO 30 J = 1,L 
JJ(J)=J 

0(I.J)=SQRTCCXCI,n-YCJ,n)**2+ (X(I,2)-YCJ,2))-**2) 

30 CONTINUE 
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I7.i59 



\%KITE C6, 1 HD I , (D(I ,K) , K = 1 ,L ) 

20 CUNTINUE 

DO 21 I =1 ,M ■■ 

V(KITE (6,221) I 
DO 32 J=1 ,L 

32 WRITE {b,2?2) JJ(J), D(I,J) 

31 CnNTl^JUE 

SORTING AND CI.ASSIFIC&TION BY- K-MN METHOD 
CALL SORT CJJ,D,M,L' ) ' 

ICO FORI-^AT (3112 ) 

1C5 FORf^AT C 2F10.2 ) 

106 FQRf-^AY ( 2Fri.2 ) 

110 rOR'*’'‘AT (IHl ,//,“** ^unknown SAMPLE SET 

1 * -A * * ■* *• * * -^< * * *• -*r * * y: -A * * i- ■* * * * K * * * 1*.- * * *• !<c ★ A V( •> A “ , / / , 

1 "Sample no" , sx , "x c i , n " , ax , "x ( i ,2) " ,//) 

112 format C no, 2F10,2-,// )' 

1 1 3 F 0 R !*•' A T ( 1 H J , / / , '• a a a T n a I N I '’■j G SET FOR K « -N a a a a a a a a- a 
1 BEFORE AND AFTER INFECTION aaa a aaaa?(:Aaaaaaaa 
1 "TRAINING ST", ax, ''Y(J,l)'',yx,"Y(J, 2)".//) 

lia FORMAT CI7,5F9.3,7X,5F9.3, //) 

220 format C1H1,//,"a*a K-NN CLASSIFICATION aaa aaaaaaa” ^ 

1 /, »aaaaaaEUCLIDEAn distance 0(1, J) IN 2 FACTOR aa 

1 "aaTRaIMING SET(nORMAUaa»,20X, "TRAINING SET(IMFECT) 
1 A A ” , X , " S A M P N Q " , 7 X , " B 1 " , 7 X , " B a " , 6 X , " e 7 ” , b X , B 1 0 " , 

I 6 X , " 8 1 2 'M 5 " A I " , 7 X , " A 2 " , 7 X , " A 3 " , 6 X , " A 7 " . 6 X , " A 0 % 

1 // ) 

221 format (///, "aaaSAMPs" ,13, » AAAAA”, / ) 

222 FORMAT (lOX, 10(ia, F6.3 5 , // ) 

STOP 

END 
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16,^22 U8«nUTLN£ PROGRAM POR WRITING 


SUaROUTiNE PROGRAM TOR WRITING 

subroutine: write ch^n,m ) 

DOUBLE PRECISION R 
DIMENSION R(2,36) 


BEFORE INFECTION, AFTER INFECTION AND AFTER CURE 

sample no from i to 12, 

WRITE (6,130) 

00 35 I = 1,.M 

35 WRITE C6,13l) I,CRCI,L), L=l,36,3 ) 

WRITE C6,iaO) 

00 ao I = i,M 

40 WRITE C6,131) I, ( R ( I , L ) , L = 2 , 36 , 3 ) 


00 40 I = 1,M 
WRITE C6,131) 
write (6,150) 
DO 50 I = 1 ,M 
WRITE (6,131) 


I, (R(I,U, L = 2,36,3 ) 


( RCI,L5, L=3, 36, 3 


SAMPLE NO FROM 1 TO 6 ,ILL VIRUS SAMPLE 

write (6,, 16 00 

00 60 I = i,M 

60 WRITE (6,161) I, ( R(I,L), L=i,P ) 

WRITE (6,170) 

DO 70 I = 1 ,M 

70 write (6,161) I, C R(i,U, L-10,18 3 

SAMPLE NO FROM 7 TO 12, NOT ILL VIRUS SAMPLE 
WRITE (6,260) •• 

DO 62 I = I ,M 

62 WRITE (6,161) I, ( R(I,U, L=19,27) 

WRITE (6,270) 

DO 72 I = 1,M 

72 write (6,161) I, C R(I,L), L=2S,36) - 


130 FORMAT 
1 

1 


131 format 

140 FORMAT 
1 
1 
1 
1 

150 FORMAT 
1 
1 
1 
X 

160 format 


C ////,8X , "v^*p,£pORc INFECTION, 12 SAMPLES 
” PP AK , 5 X , M - 1 , 6 X , « 2- 1 " , 6X , " 3- 1 , 
6JL, "4-i " ,6X, "5-1 " , 6X,"6‘^i'’, 6X,*'7«l'*, 6X, 
"8-1", 5X.’'10~l'', 5X,‘-11-1'*, 

5X,M2-r',// ) 

(13, 12F6o2 ,// ) 

( ////,8X,”**AFTER infection, 12 SAMPLES 

**»,//, "PEAK" ,5X, " 1 "a'’ ,6X ,»2~?%6X , '•5»2" , 

6X,"A'-2" ,6X,"5-2", 6X,"6*«2", 6X,"7«2", 6X, 
"8“2", 6X,"R«2", 5X,"10”2", 5X,"ll-2", 
SX,"l2’“2",// ) 

( ////,fiX,"^^AFTER CURE , 12 SAMPLES 

,//,"PEAK",SX, "1 -3",6X,"?-3",6X,"3-3" , 


6X, "4-3" .6X , , 6X , "6. 


6X 7-3' 


"8-3", 6X,"S-3", 5X,"lO-3", 5X^"n~3", 
5X,"l2-3",// ) 

C1H1,//,8X," ***■★*■< ILL , 6 SA7APLES 



k'dZ U8R0UTINE PROGKAM FOR WRITING 


1 "PEAK", "1-3", IfeX , 

1 "2-1", 6X, "2-2", 6X,"2-3",16X,"3-1" , 6X, 

1 "3-2", 6X, "3-3",//) 

161 FORMAT C 13, 3F9.2, lOX ,3F9.2, 10X,3FR,2 ,//) 

170 format ( ////,eXt" ** 5^*-* ILL , 6 samples **",// 

■ 1 "PEAK" ,5X, "a-i %6X, "A-2" ,6X, "il-3" , 16X, 

1 "5-1", 5X, "5-2", 6X,"5-3",16X,"6-1", 6X, 

1 "6-2", 6X,»6-3",// ) 

'2vS0 FORM AT (//////, 8 X. "**N0T ILL , 6S oMPLES" , // , "PE AK" , 

1 5X,"7-r',6X."7-?»,6X,"7-3",16X."8-l",6X, 

i "8-2" , 6X , "6-3" , 16X, "9-1 " ,6X , "9-2" , 6X , "9-3 

1 // 5 

270 F0RMATC////,8:< , "**NOT ILL 6 SAMPLES**" ,//, "P£6K» , 
1 • «X,"10-l'*,5X,‘'10-2",SX,«10-3",15X,"n-l“ 

1 5X,"l l-2",5X,"ll-3'M5X,'>12-l",5X,"l2-2", 

1 5X,»12-3",// ) 

RETURN 

END 



^22 UrsNGUTlNE PHOG^Av. FOR NORMALIZATIC).'^ 

C SUBROUTINE PROGRAM FOR NORMALIZATION 

SUBROUT JNc NOR’'-* ) 

DOUBLE PRECISiON B,3UMA 
DIMENSION Aa7,563 ,B(37,36) ,SUMA(36) 
C 

C' INITIALIZE DATA POINTS 



00 10 J - 1, N 




SUMA(J)=0.0 



n 

CONTINUE 




WRITE (6,1205 




DO 40 J = 

1, N 



DO 50 I = 

1, M 



SUMAC J) = 

SUMAC J) + 

ACI,J) 

50 

continue 




WRITE (6,121) J, 

SUMA( J) 


40 

continue 




i^'RITE (6, 

122) 



no 60 I = 

1, M 



DO 70 J = 

1, N 



6CI,J) = 

( A C I , J ) * 

100 ) / SUMAC J) 

70 

CONTINUE 




write (6,123) I, 

( S ( I , K ) , 

K = 1 ) 

6 0 

CONTINUE 



120 

FORMAT ClHl ,///, 

"i-A TOTAL 

PEAK area a * * * a * » * a * a * * * * 


i*^*:«:*7!f**’**i((*****’-«Vr* ■».'*******' A*-** 

1 '»SAMP NO” ,5'/ , "TOTAL AREA'',/// ) 

121 format ( 110, 5X, PiO.n, // ) 

122 format CtHl ,//,6X, "•^★percent NORMALIZED DATA******* 

l****A***-****-ti. y**x'>f A- A*.-** -A' ir -A- ■*■*■**** ***<:•<,-", 


1 

//, 


PEAK 

II 

% 

2X, 

M ^ „ 

1 ", a X , 

*' 1 ■" 5 


• 

ax 

f 

f 2 « 

7 n 
«r 

,ax, 

1 

”2- 

1” 


" p 

~2" 

i4X 

, "2- 

3" 

Ll '' 

, u ^ , 

>1 

3 

-1 

it 

,ax 

tr 

3- 

2" » 

1 - 

4X, 

"3 

-3% 

ax 


- 1 " 

» a X , 

"a 

-2" , 


w 

A 

,"a 

• 3" 

,a 

^ % 


1 

" 5*. 

i” 

,AX, 

”5 

-2" 

,ax 

, "S'* 

3" 

, a X , 


6 

» 1 
4 

II 

,ax 

H 

1 

6- 

2", 

1 

4X, 

"6 


ax 

,// 

^ 6X 

, "7- 

1 " 

,ax. 

rt 

7 

•* C? 

II 

,ax 

tJ 

? 

7 •“ 

3". 

1 


"6 

"1”, 

a X 

,”F. 

-2” 

, a X , 

"fl 

-3" , 

ax 

M 

t 

9 

-1 " 

• a 

^ 1 


1 

•• 9 ^ 

2” 

,4X, 

"<? 

-3" 

,ax 

, "10 

-1 

%.3X 

t 

H 

1 0 

M 


5X 

V 


4 

"10 

"3 

",3X 

<( 














1 

"11 

-1 

",3X 

II 

« 

1 i- 

2% 

3X,» 

1 1 

-3 " 

3 

V 

A 

M 

1 

1 

2-1 

1 

3X 

1 

1 

"12 

-2 

" , 5 X 

w 

« 

12- 

3 ” , 

/// 

) 











123 FORMAT CI3,iaF7.2, /, 3x, lMF7,2t // ) 
RETURN 

end 
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UBROUTINE PROGRAM FOP M£AN , S ^ 0 , , FcPCEMT COEF. VAR 


C 


c 

c 


c 


SUE^ROUTINE PROGRAM FOR ME AN , S . Q. , PERCENT COEF. VAR. 
SUf^ROUTINE CALC ( H , L ^ N , NN , KL , M ,ME ANB , SOP. , P TC V ) 
DOUBLE PRECISION SUMB , SfJB , ME ANB , PTC V , SSNB , BF TR ,B 
DIMENSION BC37,3 S) ,mEANB( 37) ,SDBC37) ,PTCV(37) , 
i SUMB(370. SSQBC37) ,BFTR(37) 


INITIALIZE DATA POINTS 
DO 10 I = 1,M 
SUMP (I) =0,0 
SSQBdJrO.O 

10 continue 


y*.'RITE (6,1241 
DO 60 I 
DO 9 0 J 
SUmbCI) 
SSQBdl 
90 CONTINUE 

ME A. MB ( I ) 
BFTRCn 
SDBCn = 
PTCVvd) 
WRITE C6 


1 

L ,N,KL 

SUMBdl + Bd,Jl 

ssoocn + Bd,j) 

= SUMBCI) / NN 
SORT (MN-.-SSQBCI) 
BFTftd) / NN 
C SOBCI) 100 ) 
125) I, ME/.N6CI), 


**2 

- SUMBCI)**2 ) 

/ MEANS Cl) 
SD»C2), PTCVd) 


80 

124 FORMAT 
1 

I 

125 FORMAT 
RETURN 


CONTINUE ■ 

C///,5:^, "**R£PP0DUCIBILITY OF SELECTED PEAKS 
MX, ''PEAK NO" ,6X,"MEAN>M0X,"S.U." , 
ex, "PT-C.V. " ,//) 

C 10X,I5, 3F14.2, // 3 


END 



o o 
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16.C24 USROUTImE pro ^Am FQH’ ME A N , S . 0 , - ;RCE-‘JT COEF, VAR, 

C SUBROUTINE PR ,AM FOR mEaN,S,D ZRCENT COEF. VAR, 

SUBROUTINE CA CR ,N ,NN ,Kt. ,M,f ''.H , SOB , P T C V ) 
DOUBLE PRECIS 3UMH , SOB , ME a Ufs :V,5SQB,BFTR 

DIMENSION R(? ) ,^EANB(23 tSDBC ^TC V ( 2 ; , SUM‘3 ( 2 ) , 

1 SSOBc BFTR(2) 

initialize DAT INTS 

DO 10 I = U 
SUM8(n=0.0 
SSQBCn-o.O^ 

10 CONTINUE 

WRITE ,12-^) 

00 ao = i,M 

DO 90 s L,N,KL 

IF ( ,J) ,GT, 18.00 0 ro .30 

IF ( ,GT. 9.0'. GO TO 30 

' SUMB( = SUMDCn + RC' 

ssupc ‘ = ssnacn + r( **2 

GO TO 

30 NN = . J 

90 CONTI 

MEANBCil = SUMO ( I) / NN 

BFTRCn = SORT (MN*SSDBCI) - SUMB(;)**2 I 
SOB (I) = SFTRCn / NN 
PTCVCn C SDB(I) ★ 100 ) / MEANS (15 
WRITE ( 6,1255 I, MEAnBCI), SD8 C 1 5 , PTC V C 1 5 
80 CONTINUE 

12« format .( IHI ,///, ''^^^RFP.PnoUCIBILITY OF SLED, PEAKS 

1 *i*>f**))r)‘^*Tlri»:**?r*)ir*****Vr**v**tk'**wVt:Sf*:*******'**:jf)>r:*r** 

1 **”,//, 11X4''PEAK NO" ♦6X,*'mEAN%10X-,'*S,D," , 

1 8X,"PT-C.V.",//5 ' 

125 format (/,10X,l5,3Fia.3,/ 5 
RETURN 
END 



o r» o 
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16.423 UBROUTI.^F. PROGRi.M FOR CALCULATION OF PEAK RATIO 

C SUBROUTINE PROGRAM FOR CALCULATION OF PEAK PATIO 

SUBROUTINE RATO(B,N,M,RD 
DOUBLE PRECISION tt«R 
DIMENSION B(7,’36) ,RC7,36) 

C 

DO 10 K = 1,M 
DO 20 I = 1,M 
DO 50 J = 1 ,N 
R ( I , J ) = B C I , J ) / B ( K , J ) 

30 CONTINUE 
20 CONTINUE 


sample nq FROM 1 TO 12, PEAK RATIO- 
BEFORE InFECTIQN.aFTER INFECTION AMO AFTER CURE 
WRITE C6,1S1) K 
WRITE (6,130) 
on 35 I = 1,M 

5S WR]TE (6,151) I,(R(I,L), L=l,56r3 ) 

WRITE (6,151) K 
WRITE (6,1/40) 

DO ii;0 I- = r,M 

40 WRITE (6,131) I, CHCI,L), L=2,36,3 ) 

WRITE (6,151) K 
WRITE (6,150) 

DO 50 I = 1,M 

50 WRITE (6,131) I, ( P(I,L), L=3, 56, 3 ) 

C 

C sample no FRO/^ 1 TO 6 ,PEAK RATIO 
C NOT ILL VIRUS samples 
WRITE (6,151) K 
wRITE(6,l60) 

DO 60 ! = 1 ,M 

60 WRITE C6, 161 ) I , C PCI ,L5 , L=1 ,9 ) ' ' 

WRITE (6,170) 

00 70 I = 1 ,M 

70 WRITE (6,161) I, C RCI,L), L=10,ia ) 

C 

C SAMPLE MO FR6-M 7 TO 12, PEAK RATIO 
C ILL VIRUvS Samples 

WRITE (6,151) K 
WRITE (6,260) 

DO 62 I = 1 ,M 

62 WRITE (6,161) I, C R(I,L3, L=1 9,27) 

WRITE (6,270) 

DO 72 I :i 1 ,M : 

72 WRITE (6,161) I, C R ( I , L ) , L=23 , 36 ) 

10 CONTINUE 
C 

130 FORMAT C :////, 8X, »^*eEFORE INFECTION, 12 SAMPLES 
1 *A",//,«FEAK»‘,5y,"j-i",6X,''2-r’,6X,”3-l'', 

1 6X,’»4-i‘: ,6X, "5-1 , 6X,*'6-1" , 6X ,"7-1", 6X 



«8-i « 


6X,"9-l% 5X,"10-1'', SX,”!!-!**, 



* 



423 UGROUTINE PROGRAM FOR CALCULATION OF PEAK RATIO 


150 format 

1 

4 

i 

1 

1 

151 FORMAT 


160 


1 5>:,"i2-r,// ) 

131 FORMAT CI3, 12FR.2 ,// ) 

140 FORMAT ( ////»6X, "*-*AFTER INFECTION, 12 SAMPLES 

1 * * *' , // , " PE A K , 5X , " 1 "2 " ♦ 6X , " 2-2 , 6 X , “3-2 " , 

1 &X,''4-2“,6>:,“5-2'’, 6X,“6-2“t 6X,"7-2", 6X, 

1 "8-2“, 6X,»R-2'', 5X,"lO-2", 5x,''n-2'\ 

1 5X,"l2-2",// ) 

C ////,8X, “-Ai-AFTER CURE , 12 SAMPLES 

,//, "PEAK" ,5X, “1-3'' ,6X, “2-3'’ ,6X,“ 3-3'' , 
6X,“4-3“,6X,''5-3'’, EX, "6-3", 6X,'’7-3“, 6X, 
“e-3'', 6X,''R-3“, 5X,'M0-3", 5X,-'ll-3'', 
5X,"l2-3'' ,// ) 

(IHI ,//,2x, "k = '',T2,3X,“tv* peak RATIO 

)s*>!rA'!ir*********iir*-*****A*-**"<f»f*:^ ■*•*•*■*******•*■ •',//) 

FOPMi'T ( // //,8X,“ NOT ILL , h SAMPLES 

"PEAK" ,5X," 1-1 “,6X, "1-2“ ,6X', "1.-3" , 16X , 
‘'2-1 *' ,6X. "2-2" , 6X,"2-3“,16X,"3«I'', bX, 
"3-2", s8X, "3-3" ,//T 
( 13, 3F9.2, 10X,5F9.2, 10X,3FR,2 .//T 
(//////, EX," ** NOT ILL , 6 SAMPLES 
"PEAK" ,5X, "/i-1 " ,6X, "4-2" ,6X, "4-3" , 16X, 
"5-i'',6X,"5-2", 6X,"5-3’MbX,"b-l", 6X, 
"6-2", 6X,"b-3",// T 

( ////,AX,“** ILL, b samples**" ,//, "PFAK" , 

5X, "7-1" ,bX, "7-2" ,bX, "7-3" , 1 6 X , " 8- 1 , bX , 
"3-2" ,6X, '*8-3" , 16X, "R-1 " ,6X , "9-2" ,bX,"9-5" 
// ) 

(//////, ex,"** ILL, 6 samples**" ,//,"PHAK" , 
1 4X, "10-1", 5X, "10-2", 5X, "10-3", J5X,"U-1", 

1 5 X , » n - 2 " , 5 X , “ n - 3 " , 1 5 X , » 1 2 - 1 " , 5 X " 1 2 - 2 " , 

1 SX,"l2-3",// 3 

PF.TURN 
END 


161 

170 


FOR Mi AT 
FORMAT 


1 

1 

1 

260 FORMAT 
1 
1 
1 

270 FORMAT 
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U8PW0GRAM for AUT0S.CALI'"G 

C SUBPROGRAM FOR AUTOSCALInG 

SUPROUTIME AUTS (H , N ,MEANB , SDB ♦ ASLB ) 
DOUBLE PRECISION mf anB , SOB , ASLB ,8 
DIMENSION B( 2 , 36 ) ,ASLB( 2 , 36 ) , ME AUa ( ? ) , SDB ( 2 ) 



WRITE (6,1203 
DO 93 I = 1 , M 
00 95 J 1 . N 

■ A SL8 ( I , J 3 = C B (I , J 3 -me aNB C n 3 / 

SDB(I) 

95 

CONTINUE 



WR^TE (6,1213 I , (ASLBd ,K3 ,K=1 

,N) 

93 

CONTINUE 



120 FORMAT ( IHl "A*aUTOSCALED PEAK RATIO******* 


1 ****** y.'****** ************** ************1^******’ i/ / I 


I ^ 



II 

PK 

RTO" 

» 

" 1- 

1 

n 

« 

ax 

% 

II 

1-2 

t\ 

ax 

f 

"i- 

3'S 

ax 

1 

«2« 

1 

!l 

, ax 

, »2- 

2" 

,^x 


"2 

-3 

11 


ax , 

"3 

-1 

It 

,ax 

. '=3 

-2 

1 


II 

3 

-3" 

,ax. 

'»a 

-I” 


ax 

It 

S 

n 

- 

2% 

ax 

II 

a 

“3 " 

,ax 

( 

1 

»5- 

1 

II 

,ax 

»"5- 

2" 

,ax 


II c. 

-3 

II 

n 

ax , 

"6 

-1 

It 

,ax 


-2 

1 

ax, 

II 

6 

r.3" 

,ax , 

// 

,6X 


"V 

" 1 

11 


ax. 

"7 

-2 

fi 

,ax 

,’’7 


1 


II 

8 

« 1 II 

,ax, 

"6 



ax 

II 

1 

8 

•« 

3% 

ax 

, 

9 

r» ^ 

,ax 

f 

1 

MQ«. 

2 



, "9- 

3" 

,ax 


"10- 

1 

M 

»3X 

II 

t 

,10 


2% 

3 X , 



1 " 10 - 3 S 3 X, 

1 ''ii-r',3X,"li-2",3x,''ll-3",3X,"12-l",3X, 

I « 12-2%3X.'M2"3” , /// 0 


121 FORMAT (13, 18F-.7,2,/,3X, 18F7.2, //////) 

RETURN ■ ' 

END 



C SUBPOUTINP: program po« SORTIN'G! ascending sequence 
r. AND CLASSIFICATION K-NN METHODS 
SUBROUTINE SORT CJJ,0,M»L ) 

c 

DOUBLE PRECISION 0 
DIMENSION JJCIOI, 0(26,10) 

C 

WRITE (6,220) 

DO SO K = 1 ,M 

LI = L - 1 

DC) ^0 I = 1, LI 

IPl = I+l 

DO 50 J P IPl , L 

IF ( D(K,I) ,LE, 0(K,J) ) GO TO 30 
TEMP=0(K,I) 

D(K , I)r:D(K,J) 

0(N,J)=TEMP 

TEMPS JJ (I) 

jjcn-jjcj) . . 

jj(j)-tf;mp 
50 CONTINUE 
^0 CONTINUE 

WRITE (6,221 ) K 
DO ns J=1,L 

«5 write (6,222) JJCJ), D(K,J) 

DO 46 KN =1,5,2 
WRITE (6,225) KN. 

IF ( KN ,EQ. 1 ) LL=i 

IF ( KN ,EQ* 5 ) LL=3 . 

IF ( KN .£Q, 5 ) LL=5 

KNAsO, . . ■ 

KNB=0* ^ ' 

DO 47 J = 1 ,LL 
IF (JJ(J) ,LE. 5 ) GO TO 10 

IF (JJ(J) ,G£. 6 ) GO TO 20 

10 KNB=KNB + 1 
GO to ni 
20KNA=KNA+1 
47 CONTINUE 

IF ( KN6 ,GT, KNA ) GO TO I 1 
IF ( KN8 ,LT. KNA ) GO TO 1 2 
llwRITE (6,224) 

GO TO 46 
12 'ariTE (6,225) 

46 CONTINUE 

DO 46 I = 1 ,L 
4fi JJ( I) = 1 
50 GONTINUF 


220 FfiPf'-iiT ClHl, //,"** CLiSSTF ICATIQ'^J, 8G’^TI^‘G +***» 

i ■<t***H-*1^+**V+*'»-**-******-<r****A*-j,*^**.i^*^y^^*^^l1 ///) 

221 FOPJ-AT ( ///, , 13 , " A *•**♦>«, / ) ’ 

222 FOF'^AT nox, lOCiO, Ff\,3 ), // ) 

224 FOF'^'A T C 5X , ” CLAS3IFifc’D AS ' Ft'F OFE. I iuFFc T I ON ' ) 

225 FQPf^AT CSX , ” CLASSIFIED AS *rt »AFTEF INFECTION'".//! 

RETURN ’ 
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in. results and discussion 

A. Transevaporator Sampling System 

The modifications of the transevaporator technique described here, 
result in better overall reproducibility and eliminates the problem 
of the precolumn blocking during desorption due to the presence of excess 
v/ater vapor. 

The major modification is that a "headspace" trapping is performed 
before the solvent extraction step. Minor modifications are the 
substitution of glass beads for glass wool, the removal of the Porasil 
E tube and the warming of the helium during the evaporation of excess 
2-chloropropane when the transevaporator is used in the extraction mode. 
Glass beads have been found to be a more efficient trapping material 
than glass v;ool. No differences have been found between silanized and 
non-silanized beads. Heating the helium carrier reduces the evaporation 
time for 2-chloropropane from about 20 minutes to 7 minutes. n-Pentane 
was substituted for 2-chloropropane in a few runs.- This resulted in a 
lower extraction of volatiles without significantly changing the shape 
of the profiles. ' 

Due to the small volumes involved, the transevaporator technique 
is faster than a macro headspace sampling using a Tenax trap or a micro 
extraction technique. Thus a complete transevaporator run with 
collection on both a Tenax and a glass bead tube takes about 30 minutes 
which is about the same time as required for a normal headspace sampling 
of a large sample on Tenax. 
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Figure 4 shov/s a typical chromatogram of a ten standard component 
mixture on 5 % VJitconol LA-23 stainless steel capillary column gas 
chromatographic system (2.6 ppb each component in distilled v/ater). 

This mixture was used for testing the efficiency of the transevaporator 
sampling procedure. 

For the technique of profile analysis to be useful, it is essential 
that replicate runs of the same sample yield near identical chromatograms. 
The reproducibility of the entire transevaporator/desorption/chromatography/ 
integration technique was investigated by performing replicate runs of a 
100 pi aqueous sample consisting of 1.3 ppb of each of the following 
components: n-prdpanol , n-butanol , 2-heptanone, benzaldehyde and 2- 

decanone. The transevaporator was operated in the 2-chloropropane 
extraction/glass bead collection tube mode. Table I lists the standard 
deviation and relative standard deviation for both retention time and 
peak area. The standard deviation for retention time varies between 5 
and 10 seconds and appears to be independent of retention time. Thus 
retention time is an adequate parameter for defining the identity (in 
the context of profile analysis) of all but the earliest peaks. 

In the case of peak areas, it is more meaningful to compare relative 
standard deviations, which, for the compounds considered, vary between 
4.5% and 9.8% in a random manner. These values are acceptable for 
compounds in the sub ppm range and may be considered surprisingly good, 
if the degree of sample manipulation is considered. Our prime interest 
is in the profile analysis of biological fluids where components in 
"Normal" samples are likely to vary over a range of 20% to 50%. Thus 



FIGURE 4. CHROMATOGRAM OF A TEN STANDARD COMPONENT MIXTURE 
ON 5 % WITCONOL STAINLESS STEEL CAPILLARY COLUMN 
(TOO m X 0.25 mm, I .D.) 


1 . 

Acetone 

6. 

Cyclohexanone 

2. 

Ethanol 

7. 

2-Methyl -2- 
heptene-6-one 

3. 

n-Propanol 

8. 

Benzaldehyde 

4. 

n-Butanol 

9. 

2-Decanone 

5. 

2-Heptanone 

10. 

2-Undecanone 
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TABLE I. STANDARD DEVIATIONS OF RETENTION TIME AND PEAK AREA 




Retention Time 


Peak Area 

Sample 

Mean Value 
(mtn) 

Standard 

Deviation 

Relative Standard 
Deviation (%) 

Mean Value 
(counts) 

Standard 

Deviation 

Relative Standard 
Deviation {%) 

n-Propanol 

26.54 

0.11 

0.41 

9401 

916 

9.7 

n-ButanoT 

36.46 

• 0.16 

0.44 

37320 

3228 

3.7 

2-Heptanone 

41.81 

0.08 

0.19 

65868 

5559 

8.4 

Benzaldehyde 

61.56 

0.10 

0.16 

83330 

3741 

4.5 

2-Decanone 

67.40 

0.08 

0.12 

63810 

6255 

9.3 
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changes in concentration of less than 100% are not considered significant 
and a relative standard deviation of 10% is more than adequate for the 
combined sampling and chromatographic technique. No increase in 
precision is obtained by considering either relative retention time or 
relative peak area. The reproducibility of the the technique is further 
illustrated in Figure 5 which shows three replicate runs of 60 yl of the 
same normal serum using the transevaporator in the 2-chloropropane/glass 
bead mode. The reproducibility is excellent, and with the exception of 
a few areas, the chromatograms are visually identical. The reproducibility 
is similar when the transevaporator is operated in the Tenax mode. 

The linearity of overall response is good when the transevaporator 
is operated in the 2-chloropropane/glass bead mode as is shown in 
Figure 6 for a variety of aldehydes, alcohols and ketones. However the 
actual recovery (comparing aqueous sample introduced into Porasil E 
tube to equivalent amount of dry sample directly introduced to glass 
bead tube) varies from about 20% to 80% for the five compounds illustrated 
in Figure 7. The recovery varies slightly over a concentration range of 
0.3 ppb to 5 ppb but is remarkably consistent, considering that the 
actual weight of each component is in the nanogram range. Because of 
this consistency, the non-identical recovery for different compounds is 
of little consequence to profile analysis which is a comparative rather 
than an absolute technique. 

The chromatogram obtained from a given sample when using the 
transevaporator sampling technique depends on a number of factors such 
as the nature and volume of the solvent used for extraction, the volume 


FIGURE 5. THREE REPLICATE CHROMATOGRAMS OF VOLATILE 


COMPOUNDS IN 60 pi OF NORMAL SERUM. 
TRANSEVAPORATOR USED IN THE GLASS BEAD MODE 




TIME , MINUTES 






FIGURE 6. PLOT OF PEAK AREA VS CONCENTRATION OF 
VARIOUS COMPOUNDS DISSOLVED IN TOO yl 
OF AQUEOUS SAMPLE 
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FIGURE 7. PLOT SHOWING PERCENTAGE RECOVERY VS CONCENTRATION 
FOR FIVE COMPOUNDS DISSOLVED IN TOO yl OF 
AQUEOUS SAMPLE 



RECOVERY , 
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of adsorbents and the temperature and volume of gases flovnng through 
the apparatus during the trapping procedure. Thus .a standard set of 
conditions must be chosen for the complete procedure when comparing 
profiles. A blank should be run periodically to check for contamination 
of the system. 

The transevaporator technique has been applied to a large number of 
diverse samples in our laboratory. It is especially attractive for 
biological samples due to the small sample size requirements. 

Chromatograms of normal serum volatiles (Figure 5) collected using 

the transevaporator in the glass bead mode were referred to earlier in 

connection v;ith reproducibility studies. Figure 8 shows the corresponding 

chromatogram for a serum sample collected from a diabetic patient. 

About 200 peaks are present in each of the chromatograms. The heavy ends 

of the chromatograms are visually different. In addition there are 

several quantitative differences between the two profiles that are 

evident from inspection of an integrator print out. Figure 9 shows the 

chromatograms obtained from a diabetic serum using the transevaporator 

in the Tenax mode. In these profiles the relative size of the 

arrowed peaks is reversed in normal and diabetic serum. It is 

interesting that Figure 9 shows a very similar chromatogram to that of 
(38l 

Liebich ^ a^. ' '' The latter workers used 5 ml of serum. 

Differences are also found between the profiles (Figure 14) of the 
sera of influenza patients and the sera of controls. These differences 
are manifested by a greater overall concentration of volatiles in the 
infected sera and by ratios of peak areas that are characteristic of 
either the infected or non-infected sera. 



FIGURE 8. 


A. CHROMATOGRAM OF VOLATILE COMPOUNDS IN 60 yl 
SERUM COLLECTED FROM A DIABETIC PATIENT. 
TRANSEVAPORATOR USED IN THE GLASS BEAD MODE.. 


B. SYSTEM BLANK 
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FIGURE 9. A. CHROMATOGRAM OF VOLATILE COMPOUNDS IN 60 ]i1 
SERUM COLLECTED FROM A DIABETIC PATIENT. 
TRANSEVAPORATOR USED IN THE TENAX MODE. 

B. SYSTEM BLANK 
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The transevaporator technique is well suited to other body fluids 
apart from serum. Figure 10 shov/s the profile of 200 pi of saliva using 
the transevaporator in the 2-chloropropane extraction mode. The late 
peaks are small but their size could be increased four fold by increasing 
sensitivity without introducing noticeable noise. This profile, v/hile 
not having as many volatiles as a serum profile, is of interest to 
biomedical studies due to the ease of obtaining samples. A profile 
(Figure 11) was also determined for 50 pi of breast milk which, due to 
its viscosity, was first diluted with an equal volume of distilled 
water. The profile is complex and the concentration of volatiles is 
larger than in serum. It is estimated that an adequate profile may be 
obtained with as little as 5 pi of breast milk. The profile has 
virtually no resemblance to commercial cow milk which has a relatively 
simple profile. The transevaporator technique is not suitable for urine 
due to the low concentration of volatiles. 

Figures 12A and 12B show the different profiles obtained from 
samples of cognac using the transevaporator in the Tenax mode and also 
in the 2-chloropane/glass bead mode. Other samples that have been run by 
this technique include wine, cola, pressed coffee oil and brewed coffee 
(Figure 13), all of which give complex profiles. Differences between 
Robusta and Arabica coffee were found for the latter two sample types. 

B . Volatile Compounds in Virus Infected Serum 

A typical chromatogram obtained from a pooled virus infected serum 
sample is shown in Figure 14. Peaks identified in the chromatogram by 



FIGURE 10. CHROrWOGRAM OF VOLATILE COMPOUNDS IN 200 yl 


OF SALIVA. TRANSEVAPORATOR USED IN THE 
GLASS BEAD MODE. 










FIGURE n. CHROmTOGRAM OF VOLATILE COMPOUNDS IN 50 yl 
OF BREAST MILK. TRANSEVAPORATOR USED IN THE 


GLASS BEAD MODE. 







FIGURE 12. CHROMATOGRAM OF VOLATILE COMPOUNDS IN 50 yl 

OF COGNAC. 

A: TRANSEVAPORATOR USED IN THE TENAX MODE 
B: TRANSEVAPORATOR USED IN THE GLASS BEAD MODE 



tC*99 



TIME , MINUTES 






FIGURE 13. CHROMATOGRAM OF VOLATILE COMPOUNDS IN 100 yl 
OF HOT COFFEE. TRANS EVAPORATOR USED IN THE 


GLASS BEAD MODE. 
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GC-MS are sunmarized in Table II. The peak numbers in Table II correspond 
to those marked in Figure 14. We have scanned about eighty peaks in the 
chromatogram, but only about twenty four peaks could be identified. In 
the chromatogram (Figure 14) peak number 17 appears as a well defined 
single peak, but its mass spectrum, in Figure 15, indicates combination 
of two peaks; cyclohexanone and 5-methyl -3-heptanone. Five of these 
substances (peaks 5, 8, 10, 13, 28) have been identified previously in 
normal serum. 

C . Profile Analysis of Virus Infected Serum 

A total of thirty-six virus infected and normal serum samples as 
v/ell as six pooled serum samples were studied. The processing of a 70 yl 
serum sample including sampling by the transevaporator technique, 
separation by gas chromatography and data input and computation can be 
completed within a 2.5 hr period. 

The profile of volatile substances in serum is complex with more 
than 150 peaks in the chromatogram. To simplify data handling, a visual 
comparison was made of all chromatograms, peaks due to background and 
stripping solvent were ignored and 37 peaks which appeared consistently 
in all chromatograms were selected as the data base. 

It was known that the coefficient of variation for the excretion of 
the major and consistently excreted metabolites was large. Dietary 
alteration produced relatively small changes in patterns but large 
individual within-subject variations were observed. Six 
replicate analyses of a pooled serum sample was used to establish the 



FIGURE 14. PROFILE OF ORGANIC VOLATILES FROM A POOLED 
VIRUS INFECTED SERUM BY CAPILLARY COLUMN 


GAS CHROMATOGRAPHY 
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TABLE II. SUBSTANCES IDENTIFIED BY GAS CHROiMATOGRAPHY/MASS 
SPECTROMETRY IN VIRUS INFECTED SERUM. COMPONENTS 
INDICATED BY NUMBERS IN FIGURE 14. 


Peak Molecular 


No. 

Compound 

m/e 

Weight 

1 

2-methyl -1 -hexene 

57,41 ,29,40,45,58,56,27 

98 

2 

di me thyl cyclopentane 

56,41,55,70,42,69,39,57 

98 

5 

2-propanol 

45,43,59,31,27,29,44,46 

60 

7 

2-hexanone 

43,58,51,41,29,100,85,71 

100 

8 

n-hexanal , 

44,41 ,43,29,45,55,72,82 

100 

10 

n-butanol 

56,31,43,41,26,27,42,57 

74 

12 

2-heptanona 

43,53,71.59,41,114,99,85 

114 

13 

heptanal 

44,70,43,55,57,42,29,41 

114 

14 

4-heptanone 

71.43,41,27,39,42,114,29 

114 

15 

4-octanone 

43,71,57,58,85,41,86,128 

128 

16 

6-methyl -2-heptanone 

43,58,71,95,85,110,70,59 

98 

17 

cyclohexanone and 

55,43,98,27,42,69,70,4V 

98 


5-methyl -3-heptanone 

43,57,71,72,99,29.128,41 

128 

19 

2-octanone 

43,58,59,71 ,85,128,123,109 

128 

20 

n-octanal 

43,44,41,56,57,84.29,55 

128 

22 

6-methyl -5-heptene- 
2-one 

43,44,108,55,69,58,67,111 

126 

23 

5-nonanone 

51,85.58,41 ,142,43,42,100 

142 

25 

2- ethyl hexanal 

57, 72. 4] ,43,29, 56, 71, 55 

128 
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TABLE II. (continued) 


Peak 

No. 

Compound 

m/e 

Mol ecul ar 
Weight 

27 

2-octenol 

57,56,44,43,41,29,70,82 

128 

28 

benzaldehyde 

106,105,77,41,42,40,78,43 

106 

29 

2-ethyl -1-hexanol 

57,55,41,43,46,70,83,69 

130 

30 

2-decanone 

58,43,59,71 ,29,57,41 ,55 

156 

32 

: 0 -to 1 ualdehyde 

120,119,91 ,65,82,121 ,51 ,59 

120 

33 

acetophenone 

105,77,120,51 ,43,78,106,50 

120 

34 

tri methyl -2- 
cyclohexanone 

96,138,81 ,109,67,95,97,123 

138 




FIGURE 15. MASS SPECTRA OF PEAK NUMBER 17 IN VIRUS 
INFECTED SERUM CHROMATOGRAM (FIG. 14), 
STANDARD CYCLOHEXANONE AND STANDARD 
5-METHYL-3-HEPTANONE. 
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magnitude of system variation within the data base (i.e. the effect of 
experimental variables on reproducibility). All peaks were normalized 
and the percentage relative standard deviation calculated. Reproduci- 
bility depended very much on the compound itself and consequently a 
wide variation was found between individual peaks reflecting their 
different chemical properties. The smallest variation was 6.5% (peak 20), 
a median value 18.3% (peak 5) and the largest 51.1% (peak 21) relative 
standard deviation. 

Another source of variance is the variability between individuals 
in any one group. This was calculated for each group in turn for the 
37 peaks. All peaks were normalized and their variance calculated. For 
example, the 12 serum samples of Group I (normal serum) had a lowest 
value of 17.1% (peak 1), a median value 32.9% (peak '32) and a highest 
value of 87.5% (peak 4) relative standard deviation. Results for 
Group (III) were similar but Group (II) showed far less variation. 

A second data set was selected based on the magnitude of the 
difference between the averaged normalized peak areas for the different 
groups. The criteria for selection was that the relative standard 
deviation of the selected peaks must be greater than the variation 
betv/een individuals in any one group. Seven peaks v/ere found to meet 
this requirement (peaks 16, 19, 20, 25, 29, 31, 34). 

The normalized peak areas in the second data set can show both 
negative and positive variation between groups. Under these conditions 
a ratio of tv;o peaks should prove more sensitive to inter-group differences. 
Also," if there is an interaction between peaks in each group then peak 
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ratios will be more reliable for classification purposes. The seven 
normalized peak areas were arranged in ascending or.der and all possible 
peak ratios calculated for Groups I and II. The two peak ratios 16/25 
and 29/34 were found to be most suitable for sample identification. 

The reproducibility of retention time and normalized peak areas 
for the four selected diagnostic peaks in the pooled serum sample (six 
replicate analysis) is shown in Table III. Retention times can be 
reproduced very accurately in the analytical system and this forms an 
adequate parameter for peak identification. The relative standard 
deviation of the normalized peak areas between individuals for the four 
selected peaks is’given in Table IV for Groups I through III. The 
variation in peak area for infected serum (Group II) is much less than 
normal serum (Group I) and can be more correctly defined. Peak 25 in 
Group I shows a greater variation than the others due to the inclusion 
of one extraordinary large peak in the data set. 

1. Differentiation Between Normal Serum and Virus Infected Serum 

To differentiate be tv/een normal serum (Group I) and infected serum 
(Group II) the two peak ratios 16/25 and 29/34 are calculated in Table 
V. Visual inspection shov/s that generally the peak ratio 16/25 usually 
decreased upon virus infection (except samples 7, 8, 9) and that the 
peak ratio 29/34 increased by virus infection (except samples 2, 3, 5). 
However, the use of either peak ratio does not in itself provide a 
sufficient classification between the two groups. The autoscaled data 
set of Groups I through III is plotted in two peak ratio dimensions in 
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TABLE III. REPRODUCIBILITY OF NORMALIZED PEAK AREA AND RETENTION TIME 
WITHIN A POOLED SERUM SAMPLE (6 REPLICATIONS) 


Peak 

Number 

% Normalized Peak Area 

Retention time (min) 

mean 

S.D. 

% C.V.* 

mean 

S.D. 

% C.V.* 

16 

5.36 

1.44 

26.9 

53.75 

0.57 

1.0 

25 

1.34 

0.25 

18.8 

65.82 

0.63 

0.9 

29 

4.15 - 

1 .37 

31.7 

76,51 

0.62 

0.8 

34 

1.42 

0.29 

21.0 

89.77 

0.80 

0.9 


percent relative standard deviation 
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TABLE IV. VARIATION OF NORMALIZED PEAK AREAS FOR THE SELECTED 
FOUR PEAKS IN ALL SERUM SAMPLES 


Peak 

Number 


% Relative Standard Deviation 

Before Infection 
(Group I) 

One Day After 
Infection 
(Group II) 

14-21 Days After 
Infection 
(Group III) 

16 

48.9 

32.6 

44.2 

25 

89.9 

20.4 

47.4 

29 

36.3 

32.7 

30.5 

34 

53.8 

30.3 

62.9 



TABLE V. SELECTED TWO PEAK RATIO DATA AT DIFFERENT VIRUS INFECTED CONDITIONS 


Sample 

Number 


Peak Ratio (16/25) 



Peak Ratio (29/34) 


Before 

Infection 

One Day After 
Infection 

14-21 Days After 
Infection 

Before 

Infection 

' One Day After 
Infection 

14-21 Days After 
Infection 

Serum Samples with 

Clinical Symptoms After Exposure 


- 


1 

17.73 

4.18 

8.73 

3.53 

7,30 

5.54 

2 

T7.93 

8.90 

9.65 

3.85 

2.56 

3,11 

3 

24.28 ■ 

4.77 

3.71 

44.00 

3.24 

6.17 

4 

3.32 

2.80 

3.67 

3.57 

4.65 

1.88 

5 

15.78 

2.6t 

5.77 

4.15 

3.70 

3.55 

0 

17.48 

9.51 

12.07 

0.69 

. 2.43 

1.82 

serum Samples with 

no Clinical Symptoms 

After Exposure 

i 



7 

3.38 

9,03 

8.97 

0.99 

4.13 

1.35 

8 

6.58 

5.05 

31.51 

0.66 

4.40 

5.21 

9 

1.35 

1.44 

7.56 

6.13 

7,23 

23.22 

10 

9.24 

4.95 

2.75 

1.94 

11.89, 

1.75 

11 

33.69 

5.59 

4.59 . 

1.54 

2.52 

4.00 

12 

17.59 

8.24 

23.15 

1.99 

8.43 

2.22 
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Figure 16. It can be seen that two separate clusters are formed for 
Groups I and II and thus the two peak ratios selected are adequate 
for the identification of normal and infected serum. 

In order to test the predictability of the proposed method the K- 
nearest neighbor approacii was used. A randomly chosen training set of 
ten serum samples (5 normal and 5 Infected) v/as used to assess the 
predictive accuracy of the two peak ratio technique using the remaining 
14 samples as unknowns. Figure 17 shows a computer print ous of 
Euclidean distance in two factor space (16/25 peak ratio and 29/34 peak 
ratio). Five normal sera training set is designated as 1 5 (B^ 8^2) 

and the other five infected sera training set as 6 10 (A-j Ag). In 

a given sample, the Euclidean distance between sample point and ten 
known points of the training set were rearranged by a sorting process. 
Figure 18 is a final print out for K-NN classification "before 
infection (normal)" and "after infection". This can be done simply by 
reading the first columns in Figure 18, according to already fixed K 
values (k = 1, 3, 5 ...etc.). ■ ’ 

Predictive ability v/as calculated as the percentage of the 14 . 
unknowns correctly classified. The three nearest neighbor (3-NN) 
computation assigned 85.7% of the samples correctly (1-NN, 71.4%; 5-NN, 
78.5%). Clearly the proposed method has excellent selectivity for the 
differentiation of normal and virus infected serum. 

After 14-21 days (Group III) complete recovery from infection is 
observed clinically. However, the two peak ratio classification places 
approximately half the serum samples in the virus infected category. 



FIGURE 16. AUTOSCALED TWO PEAK RATIO DIAGRAM. 

I. Normal Serum before infection; 

II. Infected sera (24 hr); 

0 , Recovered sera (14-21 days after 
infection). 



AUTOSCALED PEAK RATIO 





FIGURE 17. COMPUTER PRINT OUT FOR EUCLIDEAN DISTANCE OF 
26 SAMPLES USING 10 KNOWN SAMPLES AS A 
TRAINING SET IN TWO PEAK RATIO FACTOR SPACE. 
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FIGURE 18. FINAL COMPUTER PRINT OUT OF K-NN METHOD TO 


DIFFERENTIATE BEFORE AND AFTER INFECTION GROUPS 
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The implications of this finding can only be speculated on at the moment, 
especially as normal samples are correctly classified by this technique. 

2. Differentiation Between Normal Serum with Clinical Symptoms 
After Exposure and Normal Serum with no Clinical Symptoms 
After Exposure 

The general usefulness of this profiling technique would be enhanced 
if it could be used to predict susceptibility to virus infection. The 
normalized two peak ratio data for Group I (normal serum) samples is 
shown in Figure 19. Samples II through 16 were normal serum samples 
from volunteers who upon exposure to virus infection developed clinical 
symptoms of the disease and samples 17 through 112 did not develop 
clinical symptoms upon infection. Again two well defined clusters (with 
two exceptions 16, 19) are formed and demonstrate the possible use of 
the method for the diagnosis of virus susceptibility. However, the 
twelve samples available are too small a data base to provide a training 
set and sufficient unknowns to test the predictability of the method. A 
much larger sampling program will be required to assess the accuracy of 
this method for the determination of susceptibility to virus infection. 



FIGURE 19. TWO PEAK RATIO DIAGRAM. SAMPLE II THROUGH 

16 ARE NORMAL SAMPLES FROM VOLUNTEERS WHO 
CONTRACTED VIRUS INFECTION AFTER EXPOSURE. 

17 THROUGH 112, AS ABOVE EXCEPT NO CLINICAL 
SYMPTOMS DEVELOPED AFTER EXPOSURE. . 

0 , Recovered sera (14-21 days after infection) 
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IV. CONCLUSIONS 

Complex profiles may be obtained from aqueous samples in the volume 
range from about 5 pi (for very concentrated samples) to about 200 yl 
by the transevaporator technique for a wide variety of biological and 
other sample types. Both a "headspace" and an extraction profile may be 
obtained from the same sample. These methods are complementary as 
each emphasizes different parts of the profile. The precision of the 
method is good considering the small sample and the trace quantities of 
volatiles' present. 

The transevaporator sampling technique is shown to be useful for 
the volatile profile analysis of 70 yl of influenza virus infected sera. 

A two peak ratio method has been developed for the characterization of 
normal and virus infected sera with a percentage predictability of 85.7% 
of correctly classified unknowns. A similar two peak ratio method is 
indicative as a possible means for the assessment of virus susceptibility. 

The technique may be used to show differences between normal and 
diabetic sera, and betv^een Arabica and, Robusta coffee for both pressed 
coffee oil and brewed coffee. 



72 


V. -BIBLIOGRAPHY 

1. Biqgers, R. E., Hilton, J. J., and Gianturco, M. A., J. Chromatoqr. 
Sci . , 7, 453 (1969). 

2. Hawkes, S. J.» and VJheaton, M., J. Food Sci . , 32 , 629 (1967). 

3. Moskowitz, H. R. , J. Agr. Food Chem ., 22 , 740 (1974). 

4. Blumenthal, M. M. , Trout. J. R. , and Chang, S. S., J. Am. Oil Chem. 

Soc ., SZ, 496 (1976). • , 

5. Dupuy, H. P., Rayner, E. T., and Wadsworth, J. I., J. Am. Oil Chem . 
Soc ., S3, 628 (1976). 

6. Bullard, R. W., and Holguim, G. , 6. Agr. Food Chem ., 25 , 99 (1977). 

7. Parliment, T. H., and Scarpellino, R., J. Agr. Food Chem . . 25 , 97 

(1977). \ ’ 

8. Nordby, H. E., and Nagy, S., J. Agr. Food Chem ., 25 , 224 (1977). 

9. Palamand, S. R., and Hardwick, W. A., Tech. Quart ., 6_, 117 (1969). 

10. Dirink, P., Schreyen, L., and Schamp, N., J. Agr. Food Chem ., 25 , 

759 (1977). 

11. Grob, K., and Grob, G. , J. Chroma togr. Sci ., 8, 635 (1970). 

12. Nurok, D., Anderson, J., and Zlatkis, A., Chroma tographi a , 11 , 188 
(1978). 

13. Bertsch, W., Chang, R. C. , and Zlatkis, A., J. Chromatoqr. Sci ., 12 , 
175 (1974). 

14. 0veri;on, E. B,, Bracken, J., and Laster, 0. L. , J. Chromatoqr. Sci . , 
1^, 169 (1977), 

15. Lee, H. L., Bartle, K. D., and Novotny, M. V., Anal. Chem. , 47, 540 
(1975). 

16. Kawahare, F. K. , 0. Chromatoqr. Sci ., 10 , 629 (1972), 

17. Bentz, A. P., Anal. Chem ., 48, 454A (1976). 

Flanigan, G. A., and Frame, G. M., R/D , 28(9) , 28 (1977). 


18 . 



73 


19. Zlatkis, A., Bertsch, W. , Lichtenstein, H. A., Tishbee, A., Shunbo, 

F. , Liebich, H. M., Coscia, A. M. , and Fleischer, N., Anal. Chem., 

45, 763 (1973). 

20. Zlatkis, A., Lichtenstein, H. A., and Tishbee, A., Chroma tographi a, 

6, 67 (1973). 

21. Liebich, H. M., Al-Babbili, 0., Zlatkis, A., and Kim, K. , Clin. 

Chem ., 1294 (1975). 

22. Malya, P. A. G., Wright, J. R. , and Nes, W. R., J. Chromatogr. Sci . , 

9, 700 (1971 ). 

23. Krotoszynski, B. , Gabriel, G., O'Neill, H., and Claudio, M. P. A., 

J. Chromatogr. Sci ., 15 , 239 (1977). 

24. Teranishi, R., Mon, T. R., Robinson, A. B., Gray, P., Pauling L., 
Anal. Chem. , 44 , 18 (1972). 

25. Stoner, E., Cowburn, D. , Craig, L. C., Anal. Chem ., 344 (1975). 

26. Robinson, A. B., Pauling L., Clin. Chem ., 20 , 962 (1974). 

27. Dirren, H., Robinson, A. B., and Pauling, L., Clin. Chem., 21, 1970 

(1975). ~ 

28. Robinson, A. B., Dirren, H., Sheets, A., Miguel, J., and Lundgren, 

P. R. , Exp. Gerontol ., 11 , 11 (1976). 

29. Anbar, M. , Dyer, R. L., Scolnick, M. E., Clin Chem . ,22, 1503 (1976). 

30. Chalmers, R. A., Mealy, M. J. R., Lawson, A. M. , Hart, J. I., and 

Watts, R. W. E., Clin. Chem ., 22, 1292 (1976). • ' 

31. Malcolm, R. D., and Leonards, R. , Clin. Chem. , 22 , 623 (1976). 

32. Kitagawa, T. , Smith, B. A., and Brown, E. S., Clin. Chem., 21, 735 

(1975). - 

33. Eldjarn, L., Jellum, E., and Stokke, 0., Cl in. Chem. , 21 , 63 (1975),. 

34. Bultitude, F. W. , and Newham, S. J., Clin. Chem ., 21 , 1329 (1975). 

35. Goodman, S. I., Holland, P., Stokke, 0., Flatmark, A., and Jellum, 

E. , J. Chromatogr. , 142 , 497 (1977). 

36. Jellum, E. , J. Chromatogr ., 143 , 423 (1977). 

37. Liebich, H. M. , and Al-Babbil i , 0. , J. Chromatogr ., 112 , 539 (1975). 



74 


38. Liebich, H. M. , and Wo‘11 , J., J. Chromatogr ., 142 , 505 (1977). 

39. Owen, V. M. T. , Ho, F. K. , Mazzuchin, A., Doran, T. A., Liedgren, 

S. , and Porter, C. J., Clin. Chem . , 22 , 224 (1976). 

40. Martinez, M., Nurok, D., and Zlatkis, A., Submitted to J. Chromatogr . 

41. dellum, E., Stdrseth, P., Alexander, J., Helland, P., Stokke, 0., and 
Teig, 0. Chromatogr ., 126 , 487 (1976). 

42. Lawson, A. M. , Clin. Chem. , 21 . 803 (1975). 

43. Politzer, L. R. , Dowty, B. J., and Laster, 0. L. , Clin. Chem., 22, 

1775 (1976). 

44. Brownj G. K. , Stokke, 0., and Jellum, E., J. Chromatogr., 145, 177 
(1978). 

45. Heininqer, J., Munthe, E., Pahle, J., and Jellum, E., J. Chromatogr,, 

158, 297 (1978). ■ 

46. Pfaffenberqer, G. D. , Malinak, L. R., and Horning, E. C., J. Chromatogr., 
158 , 313 (1978). 

47. Zlatkis, A., and Liebich, H. M., Clin. Chem ., 17 , 592 (1971). 

48. Matsumoto, K. E., Partridge, D. H., Robinson, A. B., Pauling, L., 

Flath, R. A., Mon, I. R., and Teranishi, R., J. Chromatogr ., 85 , 31 
(1973). 

49. Novotny, M. , McConnell, M. L., Lee, M. L., and Farlow, R., Clin. Chem ., 
20, 1105 (1974). 

50. Zlatkis, A., Lichtenstein, H. A., Tishbee, A.,' Bertsch, W. , Shunbo, 

F. , and Liebich, H. M. , J. Chromatogr. Sci . , 11 , 299 (1973). 

51. Liebich, H. M. , Al-BabiTi, 0., Huesgen, G. , and Woll, J. , Z. Anal . 

Chem'. , 279, 148 (1976). 


52. Liebich, H. M., and Huesgen, G., J. Chromatogr ., 126 , 465 (1976). 

53. Dowty, B., Carlisle, D., Laseter, J. L., and Gonzales, F. M., 

Biomed. Mass Spectrom . , 2_, 142 (1975). 

54. Dravnieks, A., and O'Donnell, A., J. Agr. Food Chem ., 19, 1049 (1971). 

55. Schultz, T. H., Flath, R. A., and Hon, T. R., J. Agr. Food Chem. , 19 , 
1060 (1971). 



75 


56. Zlatkis, A., and Andrawes, F., J. Chroma toqr ., 112 , 533 (1975). 

57. Stafford, M. , Horning, M. G., and Zlatkis, A., J. Chromatog r., 126, 

495 (1976). ~ 

58. Zlatkis, A., and Kim, K. , J. Chromatogr ., 126 , 475 (1976). 

59. Isenhour, T. L., and Jurs, P. C., Anal. Chem ., 43 , 20A (1971). 

60. Jurs, P. C., Kowalski, B. R., Isnehour, T. L., and Reilley, C. N., 

Anal. Chem ., 690 (1969). 

61. Jurs, P. C.j Kowalski, B. R., Isenhour, L., and Reilley, C. N., 
Anal. Chem. , 4^, 1387 (1970). 

62. Kowalski, B. R. , Jurs, P. C., Isenhour, T. L., and Reilley, C. N., 

Anal. Chem ., 695 (1969). 

63. Cover, T. H. , and Hart, P. E., IEEE Trans, bn Info. Theory, IT-13, 

21 (1967). 

64. Peterson, D. W. , IEEE Trans, on Info. Theory , IT-14 , 26 (1970). 

65. Kowalski, B. R. , and Bender, C. F., Anal. Chem. , 44 , 1405 (1972). 

66. Leary, J. J., Justice, J. B. , Tsuge, S., Lowry, S. R., and Isenhour, 

T. L., J. Chromatogr. Sci . , 11 , 201 (1973). 

67. Clark, M. A., and Jurs, P.- C., Anal. Chem ., 47 , 374 (1975). 

68. Lawson, A. M. , Chalmers, R. A., and Watts, R. W. E., Clin. Chem ., 

1283 (1976). 

69. Chalmers, R. A., Healey, M. J. R. , Lawson, A. M., and Watts, R. W. E. , 
Clin. Chem ., 22 , 1288 (1976). 

70. Anand, V. D., White, J. M. , and Nino, H. V., Clin Chem., 21, 595 

(1975). - 

71. Glaser, E. R., Silver, B., and Suffet, I. H., J. Chromatogr. Sci ., 
15,22(1977). 

72. Liebich, H. M., and Woll, J., J. Chromatogr ., 142 , 505 (1977). 

73. Cadogan, D. F. , and Sawyer, D. T., Anal . Chem . , 42 , 190 (1970), 

74. van Wijk, R., J. Chromatogr. Sci., 8, 418 (1970). 

75. Stenhagen, E,, Abrahamsson, S., and McLafferty, F. W., (Editors), 
Registery of Mass Spectral Data, Wiley, New York 1974. 



76 


76. Eight Peak Index of Mass Spectra, Mass Spectrometry Data Center, 

A1 derma rston, 2nd Ed., 1974. 

77. Jurs, P. C., and Isenhour, T. L., Chemical Application of Pattern 
Recognition , Wiley-Interscience, New York, p. 30, and p. 78, 1975. 


78. Jeltes, R., J. Chroma togr. Sci. , 12 , 599 (1974). 



